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2 ■ 1. INTRODUCTION 

S: 

■ Various macroscopic phenomena seen in high temperature superconduc- 
Q . tors (HTSC) under nonzero magnetic fields have inspired much interest in 
O i studying the phase diagram of vortex states of three dimensional (3D) layered 

type II superconductors. Experimentally, thermodynamic and transport phe- 
nomena have been intensively examined mainly in YBCO [1-3] and BSCCO 

■ [4]. The resistive broadening [5] in the so-called vortex liquid regime below 
<N : the crossover line Hc2{T) is common to these materials with strong fiuctua- 
i/^ '. tion effect and is well described in terms of the non-Gaussian superconducting 
O I fiuctuation theory [6]. Since, on approaching a phase transition, macroscopic 

behaviors tend to become insensitive to microscopic details of each material, 
one expects that the phase diagram itself should be qualitatively the same 
between these materials. However, recent data [7-10] in YBCO seem to have 
shown some details of transition lines in real systems with pinning disorders 
I which have been unseen in BSCCO. It is important to find how such phe- 

Q ' nomena apparently dependent on the materials are explained within a single 

O ■ GL theory. 

Theoretically, the phase diagram of vortex states has been tackled from 
^ ! two different points of view. In the literature [11-13] based on the elastic 

^ I theory for the vortex systems with pinning disorder, one first starts from low 

fields and low temperatures and hence, works in the London (phase-only) 
limit of the GL model. Consequently, a glassy solid phase, named Bragg- 
glass (BrG) phase, was proposed as the ground state in cleaner real systems. 
Then, a melting line of BrG phase is estimated in terms of some kind of Lin- 
demann criterion and is usually identified with a simultaneous destruction of 
positional and glass (superconducting) orders. Since the positional ordering 
of field-induced vortices in clean bulk systems is believed to occur through 
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a first order transition [14], however, this picture of glass ordering [11-13] 
does not permit a continuous disappearance of hnear resistance and hence, is 
incompatible with transport phenomena in fields lower than a lower critical 
point of clean YBCO samples [7-9] and in systems with a continuous glass 
transition [15] due to strong correlated (line-like) disorder. 

By contrast, the approach [16-19] from higher temperatures describes the 
vanishing of linear resistance according to a superconducting glass ordering 
proposed by Fisher et al. [20] for homogeneous (nongranular) type II super- 
conductors in nonzero fields. In layered materials, this ordering occurs when 
the glass susceptibility, which is the spatial average of correlation function 
[20] 
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expressed in terms of the pair-field (superconducting order parameter) ipj{r) 
at j-th layer, becomes divergent. In eq.(l.l), the angular bracket and the 
overbar denote, respectively, the thermal and the random averages, and d is 
the interlayer spacing. This approach can be formulated [17-19] as a natural 
extension of the nonGaussian superconducting fluctuation theory [6] to lower 
temperature at which the vortex pinning due to structural disorder is not 
negligible even in clean systems. Previously, this approach was criticized in 
a review paper [21] because it was not easy [16] for this approach to justify 
the phenomenological guess [20,21] that, in thermodynamic limit, the vortex 
solid in the pinning-free case with nonzero vortex flow resistance should be 
replaced in clean limit of real systems by a vortex glass with zero linear 
resistance. This obstacle was overcome [17,22] at least in high field case by 
combining properties of the pinning-free Abrikosov solid in 2D limit with 
the framework of vortex glass fluctuation based on eq.(l.l). Further, since 
eq.(l.l) is an expression independent of the detail of random average, this 
approach is easily extended [19] for describing continuous glass transitions 
[15] induced by correlated disorder. On the other hand, it is unclear at 
present to what extent this approach can be extended into the resulting 
glass phases. 

In this article, a theoretical development on the vortex glass transitions 
based on cq.(l.l) is reviewed. Consistently with eq.(l.l), the GL model for 
the layered system, i.e., the Lawrence-Doniach (LD) model [23], will be used 
throughout this paper: 
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with disorder terms 



Hrp^dT, I d'r 



:i.3) 



where ,^0 the in-plane coherence length, 0o the flux quantum, b > 0, T the 
mass anisotropy, and the vector indices _L and || imply the directions, respec- 
tively, parallel and perpendicular to the layer plane. Throughout this paper, 
internal gauge fluctuations except the external disturbance SA are neglected 
by focusing on the type II limit so that the applied fleld is given by curlAgxt 
where Acxt = A - 6A. In eq.(1.3), j/ = V^f (-iV + 27rAcxt/0o)^« + cc, and 
the structural disorder is described by a random potential expressing Tcq- 
variations, Uj{r), and a randomness of flux /j(r) [18,19]. We note that, in 
the phase-only hmit (|'0| const.), the second term of eq.(1.3) expresses the 
pinning of vortex cores [11,12], while its first term becomes negligible. 

This paper is organized as follows. In §2, the phase diagram of real 
systems in the case with B _L layers and, primarily with only point (uncorre- 
lated) disorder, is discussed according to recent works [17-19, 24], which are 
based on the theoretical findings that, in the pinning-free case, the vortex 
liquid region of the normal metal phase discontinuously freezes to change 
into the Abrikosov vortex solid with long-ranged positional order and quasi 
long-ranged [25] (conventional) phase coherence and that, in real systems 
with pinning disorder, a static glass ordering defined using eq.(l.l) occurs 
while the conventional phase coherence remains short-ranged because of a 
disorder-induced partial destruction [20,21] of positional long-ranged order. 
The above second statement may be subtle if the BrG phase [11-13] will oc- 
cur just below the first order freezing of the vortex liquid. In §3, the phase 
diagram in the case B || layers of real systems with point disorder is discussed 
by applying [26,27] the treatment sketched in §2 to this case and is briefiy 
compared with existing data [28,29]. In §4, a relevance of results in §2 to 
BSCCO in low fields perpendicular to the layers is discussed together with 
an issue of Hall conductivity near glass transitions. 



2. PHASE DIAGRAM IN FIELDS PERPENDICULAR TO LAYERS 

First, the model (1.2) will be rewritten within the subspace of the lowest 
Landau level (LLL) of the '?/^-fiuctuation by invoking a high field approxima- 
tion valid far from a critical region of the normal-Meissner transition at Tcq. 
Further, the 2D case will be considered [22] for a while for our convenience of 
description. Using a Laudau-gaugc and expressing as ^(r) = J2p(ppUp{r) 
in terms of the LLL eigenfunction Up, the model Ti. = Hld + Ti-rp takes the 
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form 

^ = E /^o I0pr + E(;i} iPkr + (^-k + /-k k^eo) ^i^'pk) , (2.1) 

where = /(p^, jiQ = —1 + h + T/Tco, B the magnitude of the apphed 

field, = exp(-k^/(2/i)), and 

Pk = y E exp(ipA;i//i) (/)*_fc2/2 </'p+fe2/2 (2.2) 

with hnear system size L. Since the characteristic microscopic length of a 
vortex state is the magnetic length = ^ (f)o / (2% B) , the factor k"^ of the 
random-fiux term implies that the random potential / is accompanied by the 
factor h and hence that the pinning disorder will be enhanced with increasing 
field. This trend is also valid in the phase-only model with no rt-potential 
term and will be valid in general at least in type II limit [20,21]. The Gaussian 
ensemble for the w-potential will be assumed: u{r) u{r') = AS^'^\v — r'). For 
just simplicity of our presentation, the /-potential term will be omitted for 
a while. 

First, let us consider, as typical quantities appearing even in the pinning- 
free case, the pairing entropy density (pair-field propagator) 

< > = I d?r< IV^P > (2.3) 

and the Abrikosov factor 

Pa = d'r'^W^y'j d'r< >, (2.4) 

where Ny = L^/ (27rr^) is the number of vortices. The former satisfies a 
Dyson equation 

= =J^= ^fXo + fxx{PA- A,s{T){Pa - 1)), (2.5) 

where x = k^T bh/{27i^Qdfi'^), and Acff(T) = Ad^Q/k^Th is the pinning 
strength relative to the thermal fluctuation strength. The Abrikosov factor 
[30] is expressed in the form 

= 1 + iV-i ^ ^;k( 1 - 2 X He ), (2.6) 

k 

where Vk is the fuUy-renormalized four-point vertex corresponding to the 
bare one fk and also depends on Agg. If the freezing to the sohd is of first 
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order, a precursor of the positional ordering of vortices (zero points of ip) 
will appear only in the vicinity of the transiton. Then, sufficiently above the 
transition, Vk will take a form of RPA type such as 

yiiq ^ . (2.7) 

On the other hand, one will notice by comparing eq.(2.6) with its mean field 
expression [30] that in the limit of a perfect solid takes the form 

= V^'^ix :^1) + SV^=^{1-N,Y: ^k,G ), (2.8) 

G^O 

where G's are the reciprocal lattice vectors of the vortex solid. Note that 
the first term corresponds to the x — > oo (i.e., low T) hmit of V^"^, while 
the second term 5Vk corresponds to the structure factor of a vortex state. In 
a solid-like vortex state with a long but finite positional correlation length, 
SV\i near G will take, for instance, a Gaussian form like 

where (3> 1) corresponds to the positional correlation area A^cor- One 
can verify by substituting eq.(2.9) into eq.(2.6) that the /3^-value does not 
depend remarkably on A'^cor at least at low enough T. In addition, we note 
that, at low T, eq.(2.5) reduces to the mean field result 



«|VP»sp = -(&/3a)-Vo, (2.10) 

where < >sp denotes space average, and the small O(Aeff) correction to (3a 
was neglected. 

Now, let us turn to the glass susceptibihty in 2D, which is expressed 
within LLL in the form 



XG = K\< >)-2^ I < 0,0;, > |2. (2.11) 

p,p' 

In clean limit, xg is given as a ladder-series of the irreducible vertex repre- 
sented in Fig.l, i.e., 

XG = 1 + /irr + + ^frr + , (2-12) 

where 



^'"=2^S-.(l-2xy.)^. (2.13) 
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Far above the freezing transition where Vk is dominated by the RPA term 
(2.7), we obtain ^ Aeff/2(1 — 0(x~^)). Since this expression is insensitive 
to B, one may conclude the absence [31] of 2D glass transition even at the 
mean field level. However, once the quantum 

superconducting fluctuation is taken into account, a glass transition line 
at the mean field level can exist above a (if any) first order transition line 
[32]. 

By contrast, when taking account of SV^ illustrated by eq.(2.9) which 
becomes rather remarkable in the vicinity of the (if any) 

first order line, one obtains a result suggestive of a glass transition in- 
duced [17,22] by the vortex solidification (i.e., by the first order transition). 
For instance, if substituting the expression of the perfect solid eq.(2.8) into 
eq.(2.13), one finds 

Iir,^^if3A-l)N, (2.14) 

proportional to the total number of vortices. The factor — 1 (> 0) 
implies that a spatial variation of [V'l due to the vortices is crucial in obtaining 
a glass ordering at or below the (if any) first order line. The factor is a 
consequence of the assumption of a perfect solid in the case with a small but 
finite A and must be replaced by the positional correlation area A^"cor(Acfr) 
in terms of, say, eq(2.9). Note that the origin of the large factor or A'cor 
is the vertex correction to the impurity line (the semicircles in Fig.l). If, as 
an estimation of Ncor, identifying it with the correlation area resulting from 
the collective pinning theory [33], we find Ncor — {^cs)~^- Further, since 
the mean field glass transition line Tq^{B) in clean limit will lie just below 
the freezing transition, or crossover, line Tj^'^\B) in 2D LLL, /x^ in eq.(2.13) 
may be replaced by its value at T^'^\B) IQ-'^hh {2ti ^Id)'^ . Then, since 
the /^A-value is insensitive to material parameters, the resulting is almost 
independent of material and physical parameters. At least, it does not vanish 
in clean limit (Agg 0), implying that the disorder-free theory cannot be 
used even in clean limit below the expected T^'^^B). 

Next, let us extend the above analysis to 3D case in which the glass transi- 
tion will occur more easily. For a while, the case of point disorder will be con- 
sidered in which the random potentials satisfy Uj{r)ui{r') = d~^ A^'^ S^'^\r — 
r') 6j^i and fj{r)fi{r') = d~^ A^^^ 6^'^\r — r') 6j^i. Roughly speaking, xg in 3D 
case is given by replacing A//x^ in eq.(2.13) by A^^ /j,^^"^ when, for simplicity, 
neglecting a A$ term. If, as in 2D case, A^cor is identified with the dimen- 
sionless positional correlation area perpendicular to B found in the collective 
pinning theory (in type II limit) [33], we obtain 

/i„ - A,s{Pa - 1) exp(ci A;i) (2.15) 
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along the first order line expected in LLL (see eq.(2.17) below), where ci is 
a positive constant. Since eq.(2.15) is divergent in A^e limit, it implies 
that the glass transition due to point disorder, the so-called vortex-glass (VG) 
transition, will occur more suddenly in cleaner systems with smaller Apff and 
that the first order transition in clean limit should be a glass transition 
simultaneously. Further, the exponential Aefr-dependence of A^cor suggests 
that this pinning dependence is stronger than the prcfactor Aefr in eq.(2.15). 
Since A'^cor will decrease down to a constant of order unity with increasing Agg , 
and Pa depends only weakly on A^s, it is expected that Ii„ monotonically 
decreases with increasing A^g down to a value ~ A^^h{(3A — l)/(27r/x^''^). 
Namely, the resulting VG transition line is expected to deviate from the first 
order line to Zow;er temperature with increasing Aefr(T) and, in high Aeff(T') 
limit, approach 

B^^{T) c H^2{0){^-^Y\A^^y/'. (2.16) 

Although this line [16] satisfies the LLL scaling B ~ (Tc{B) — T)^/^ as well 
as the first order line [34] in 3D and LLL 

S^(r) ^ HMO)^-^^, (2.17) 

the dependences on the (anisotropic) 3D fluctuation strength 9f — bkBTy/T/^^ 
of ByQ and B„i arc different from each other. 

Through the above findings on the glass transition at or below Bm{T), a 
picture on the phase diagram above a lower critical point Bif^ (defined below) 
of the thermal first order transition is easily obtained [17] . First, since most of 
the above expressions depend not on A^'^ but on the relative pinning strength 
Aefr(T) oc 1/T, a cooling along Bm{T) (i.e., an increase of B) will imply 
an effective enhancement of random pinning effect. Further, as mentioned 
below eq.(2.2), an inclusion of nonzero A$ additionally induces a pinning 
effect enhanced by an increase oi B [18,19]. For these reasons, the first order 
transition along B^(T) will be weakened with increasing B and will disappear 
at some upper critical point -Bucp- Then, it is clear that the vortex state 
just below the thermal first order transition and below B^^p need not have a 
positional order. Actually, as explained above, the glass transition line begins 
to deviate from Bm{T) to lower temperature with increasing Aefr(T') or B and 
approaches Byq{T) at high enough fields. Note that ByQ{T) decreases with 
reducing the pinning strength and, as suggested below eq.(2.17), decreases 
more rapidly than B^iT) with increasing Of. Hence, there is a possibility 
[17] of a wider window of the so-called vortex slush regime [35] in cleaner 
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systems as far as it is not masked by the BrG phase which may exist at lower 
fields (see Fig. 3 below). 

Of course, the vortex slush regime is a part of the vortex liquid region 
and hence, of the normal metal phase because the linear resistance is finite 
there. In the present theory, the in-plane resistivity p^x (= Pyy) in this 
regime vanishes algebraically on approacing the glass transition, as in the 
strong disordered case, but with a field-dependent smaller exponent [18]. To 
show this, let us briefly explain how to evaluate the conductivity near a VG 
transition. As accepted even through the studies [6] at higher temperatures, 
the conductivity cTjj may be separated into the quasiparticle part cr^^jj and 
the superconducting fluctuation part (Js,ij, and, deep in the liquid regime of 
clean systems, as, ij can be expressed as a sum of the pinning-free contribution 
(jp^ij and the glass fluctuation part (TG,ij so that aij ~ crn,ij + <^F,ij + <^G,ij [17]. 
The dynamics of ip-field is incorporated according to the TDGL equatioin, 
or equivalently the quantum TDGL action [17,36] 



'QLD 



/ rfV/3^(7i|a;|+i72C^)|V'.-(r,u;)|2 + £^ ±\^nU^ ^ ^l^(r)) 



(2.18) 



for the LD model, where P = 1/(A;bT), 71 > 0, and u denotes Matsubara 
frequency. As flrst found in Ref.37, the vortex flow conductivities 

CTF,.. = R;'^{-^^o) (2.19) 

and (TF,xy = l2<^F,xx/li ai'e obtained as the low T limit of the renormalized 
Aslamasov-Larkin (AL) fluctuation conductivities [38], where Rg is the re- 
sistance quantum 7:h/2e^. According to eq.(4.1) of ref.6, the corresponding 
diagonal AL conductivity parallel to B has the following low T form obeying 
the LLL scaling 

aF,..^a,,xx[—^^^) ^^sk^f^- ^^-^O) 

The same relation was derived later in ref.39 in terms of the phase-only 
model. It suggests that the behavior (2.20) of pzz is also vahd in the hquid 
regime in lower fields. 

Below, we focus on the glass fluctuation term (Tg,xx- The Feynman dia- 
grams expressing (Jg,xx are illustrated in Fig.2. In weak enough pinning case, 
the resistive vanishing can be described just by Fig.2 (a), while Fig.2 (b) 
becomes necessary in order to derive the universal VG scaling. 
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Actually, assuming the glass transition to be continuous, Fig. 2 (a) gives 

[17] 

(2.21) 

while the diagrams such as Fig. 2 (b) result in the scaling behavior [20] 
argued by Fisher et al. 

^U-iT-Tyo)^'-'^", (2.22) 

where Tvg(-B) is the VG transition line (in type II limit), and z (> 4) and u 
are, respectively, dynamical exponent and the exponent of correlation length 
Cvg(^) in VG critical region (in type II limit). aG,xx is given by a sum of 
eqs.(2.21) and (2.22). A crossover between the behaviors (2.21) and (2.22) 
occurs when [18] 



Namely, if Lcr{B) is beyond the system size, the behavior (2.21) may be seen 
like a true critical behavior in type II limit in cleaner systems, and, in a 
clean sample, the exponent of vanishing resistivity is i?-dependent and will 
increase from ^{z — 3) to ^{z — 1) with increasing B in an apparently con- 
tinuous manner [18]. In fact, a S-dependent exponent was observed in the 
vortex slush regime of a YBCO sample, and the expected universal behavior 
(2.22) seems to have been found in higher fields than -Bucp [40]. A similar 
behavior was observed previously in other experiments: the resistivity expo- 
nent in a moderately disordered sample was S-dependent and smaller than 
the expected one (2.22) [41], while a 

S-independent scaling behavior has been observed in dirtier samples [42]. 
This expectation [18] of an algebraic and 5-dependent scaling of resistance 
in the vortex slush regime is different from the thermally-activated vanishing 
in the slush regime argued by Worthington et al. [35]. 

So far, our discussion has been limited to the field range in which the 
VG transition occurs (at or) below the first order line or its extraporation 
to higher fields than -Bucp- Since, as mentioned above, the pinning disorder 
becomes effectively weaker in lower fields, it may be natural to expect the first 
order transiiton not to terminate at lower fields. However, the thermal first 
order transition should not occur any longer in fields where a glass transition 
line exists above Tm{B), and a lower critical point Sicp of the first order line 
should appear if the glass transition lies above Bm{T) in lower fields [19]. 

Actually, this situation is realized even at weak disorder in the present 
LLL approach which may be qualitatively valid far above the Hci{T). To 
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BA^s{T) V 4/i2A$ 



(2.23) 
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see this, let us first consider the case with only Zine-like disorder parallel 
to B [19,43], defined by Uj{r)ui{r') — A(^)5*^^^(r — r') and a similar one for 
the /-potential. The resulting glass transition due to such correlated defects 
II B is called in the literature the Bosc-glass (BG) transition [44]. Assuming 
the BG transition to occur above Bm{T), the 5Vk-contribution to the vertex 
correction to the impurity line, carrying A^'^ in this case, in Fig.l can be 
neglected, and the vertex correction consists only of the RPA term (a 3D 
version of eq.(2.7)). Then, Ii^j. is easily obtained in the lowest order in A*^'\ 
and a BG-line results in. When Ip^ = ^q/ \/T]l (the usual phase coherence 
length II B) 3> (i, it is expressed as 

i?BG(T) ^ ife2(0)^^-^(-/xo), (2.24) 

which is hnear in Tco — T near T^o, as observed in ref.15. Since B^{T) (2.17) 
has a vanishing curvature near T,.q, one finds that the resistivities vanish on 
cooling continuously on B-qq(T) above Bjn{T) m. B < B^^^^, where 



(AW^3 



72777^^02(0), (2.25) 



increasing with reducing the fluctuation or enhancing the pinning. Accom- 
panying this BG transition, not only the diagonal conductivities a^x and a^z 
but also the tilt modulus (the diamagnetic susceptibility to a transverse field 
± B) show a critical divergence, implying the presence of a transverse Meiss- 
ner effect in the BG phase [44] . The details of calculations of these response 
quantities near a BG transition will not be given here and can be found in 
ref.19 together with the corresponding results in the Gaussian splayed-glass 
transitions [45]. 

Interestingly, the corresponding situation with ByG(T) > B^{T) also 
occurs in the case [24] with only point disorder. Under the same assumption 
as that used above for Bbg{T), we find a ByQ{T)-lme [17] 

linear in — /xq, where Bdc{T) is the so-called decoupling crossover line [34]. In 
this case, the resulting lower critical point B^^^ 

Bii^nm{^y{^ (2.27) 
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is usually much smaller than but has similar dependences on the pinning 
and fluctuation strengths to the corresponding B^^^. We note that the above 
expression was derived in the lowest order in A^p\ An inclusion of the next 
order contribution to Ii„ tends to increase i?[^^p-value, although the resulting 
dependences on the pinning and fluctuation strengths become complicated 
[46]. Hence, cq.(2.27) should be seen as a lower limit of the expected lower 
critical point. In any case, the resistivities in B < are expected to vanish 
continuously at the second order VG transition, and the thermal first order 
transition should not occur in these low fields where Bm{T) lies in the glass 
phase. It is important to note that, as clear from the above discussion, a lower 
critical point is not due to an enhancement of pinning effect accompanying 
a lowering of the field: A lower critical point was observed in a couple of 
experiments [7-9] in tesla range of YBCO clean samples where the type II 
limit neglecting fluctuations of flux density is safely valid. It is theoretically 
difficult to expect a pinning-cnhanccmcnt accompanying a field-lowering in 
type II limit. Actually, for this reason, the appearance of a lower critical 
point was not predicted from treatments based on the vortex elasticity. 

As a test of the present explanation [17,24] on the existence of the vortex 
slush regime and of i?icp, it is interesting to compare the above results with 
the oxygen- deficiency dependence [9] of YBCO phase diagram. As systemat- 
ically examined by Nishizaki et al., both the upper and lower critical points 
of first order line tend to decrease with underdoping [9,10]. It is well known 
at present through the doping dependences of penetration depth [47] and 
heat capacity jump [48] that the thermal superconducting fiuctuation is en- 
hanced with underdoping. On the other hand, effects of point disorder due to 
oxygen deficiency also become more remarkable with underdoping. Namely, 
both 9f{Tco) and A^^^ increase with underdoping. It is natural to interpret the 
doping dependence of B^cp as being due to a Aeff(T) ~ A^y9f{T)-mcrease 
(i.e., a relative enhancement of pinning) with underdoping. However, addi- 
tional dependences on the anisotropy F and on 6f of B^^^ (2.27) imply that 

i?|^p can decrease with increasing A^s, because F remarkably increases with 
underdoping [48]. An explanation of other experimental findings on the lower 
critical points was given in ref.24. Now, we are in a position of discussing 
possible phase diagrams under B _L layers, which are described in Fig.3 [24]. 
Before explaining Fig.3, we need to mention a bit about results of the elastic 
approach. As emphasized in, for instance, ref.l2, the melting line of the BrG 
phase should be of first order in general. However, the superconducting tran- 
sition observed in B < Bicp is of second order, and it is difficult to explain 
the presence of Bi^p consistently with the argument favoring the BrG phase. 
As far as the lower critical point is unrelated to the BrG melting, there is no 
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reason why all of the thermal first order line in Bicp < B < B^cp is included 
in the BrG melting line. Through some consideration including the above 
statements, we have concluded that a generic 3D phase diagram in B ± 
layers with intermediate strengths of fluctuation and pinning will be of the 
form Fig. 3 (a), in which the BrG melting is separated from the thermal first 
order transition occuring along Bm{T), and the vortex slush regime exists 
entirely in B\cp < B < By^cp- An evidence of a BrG melting line lying much 
below Byq{T) was recently found in data of NbSe2 [49] and (K, Ba)Bi03 [50]. 
However, in HTSC with strong fluctuation, the thermal first order line may 
be pushed down to lower temperature and, in part, merge the BrG melting 
line. Then, the only possible phase diagram will be of the form Fig. 3 (b), in 
which the glass phase just below the first order transition in B\cp < B < B* 
is BrG, and the vortex slush regime exists only in B* < B < B^cp- A recent 
magnetization measurement in heavily overdoped YBCO [51] seems to have 
shown the presence of the BrG melting line in i? < B\cp approaching T^o with 
decreasing B just hke what we expect through Fig. 3 (b). The dashed curves 
in both Fig.3 (a) and (b) indicate S^(T). In passing, we note that the 
phase diagrams given in Fig.3 are valid even for real systems including weak 
line-like disorder, although in this case the glass phases have the transverse 
Meissner effect, and the vortex slush regime is much narrower [24] than in 
the case with no line disorder. 



3. PHASE DIAGRAM IN FIELDS PARALLEL TO LAYERS 

In this section, we briefiy explain results found in LLL approach to the 
phase diagram of model (1.1) in B || layers. In this field configuration, the 
field strength is usually measured by the combination [26,52] 

27rd2 

P = -—^B, (3.1) 

as far as the relation < is satisfied. Below, we primarily focus on 

the strong field region satisfying exp(— p) <^ 1 in such a layered system. In 
such high fields, the mean field transition line T(.{B) approaches a limiting 
behavior [53] 

T,(S)^T,o(l-^) (3.2) 

independent of p, and the action (2.18) written in terms of LLL modes takes 
the simple form [26,27] 

^QLD „.,o f^f'dT 



h 
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x\m,r)\' + 



b 



^o{ni - 712, n2 - na; qy^i) 

Qi,Q2,Q3 



2dL2 



x0*(Qi,r)0*(Q2,r)0(Q3,r)0(Qi + Q2 - Q3,t) , 



(3.3) 



where is the Fourier transform of the LLL fliictuation field 0(t), Qj = 
QxX + Qji), Qj == qy,j + r]^^dnj with integer rij, fiQ = {T - Tc{B))/Tco with 
eq.(3.2), and the disorder energy term Tirp was dropped for convenience of 
our presentation. The bare vertex Vq is an even function [26] of rii — ris and 
n2 — ns, implying that the partial LLL degeneracy of degree A^^ = L/d is 
measured by Uj. Further, the gauge disturbance SA necessary in deriving a 
conductivity at a uniform current was assumed to be spatially uniform. The 
corresponding action useful in deriving tilt modulus can be seen in ref.27. 

In the disorder-free case, the only true transition is argued again to be a 
first order transition at Tm,{B) between a vortex solid and a (narrow) vortex 
liquid regime below Tc{B) [26]. In low fields, Tm{B) is close to the corre- 
sponding one 

of the anisotropic 3D GL model, which is given by eq.(2.17) with F re- 
placed by 1, while it approaches, as well as Tc{B), a p-independent value in 
large p limit [26,27]: 



with eq.(3.2), where 9j: (T) = b k^T/ {271 ^^d) is the fluctuation strength in 
2D, and a constant of order unity Cm (> 0) has not been determined analyt- 
ically. The fact that Tm{B) becomes independent of p for high p values is a 
reflection of conflnement of vortices between all interlayer spacings. Because 
an increase of p in p > 1 does not delocalize the vortices out of interlayer 
spacing any longer but just compresses each vortex row along the layers, 
the spatial variation of \ip\ on the superconducting layers diminishes with 
increasing p. Namely, since the Abrikosov factor /?^, eq.(2.4), approaches 1.0 
with increasing p (> 1) irrespective of the vortex lattice structure at lower 
temperatures, the flrst order transition becomes signiflcantly weaker with in- 
creasing p. In the simulation of 3D XY model [54] where Tc{B) is always 
identical with Tgo, a similar melting line insensitive to p in p > 1 was detected 
from heat capacity data. 

However, the high p (> 1) portion of the melting transition was argued 
there not to be weakened flrst order mentioned above but to be continuous. 
This controversy may not be resolved by real experiments because, as dis- 
cussed below, there is a reason why the disorder-free melting transition in 
higher p should be easily destroyed by disorder existing in real systems [27]. 




) 



(3.4) 
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Here we merely note that of eq.(3.4) decreases with increasing the fluc- 
tuation strength 9f^\Tco) or with decreasing the anisotropy F. The latter 
dependence seems to become dominant in the doping dependence in YBCO 
according to resistivity data in ref.55. 

Below Tm{B) of the disorder-free system, a Josephson- vortex-solid, i.e., 
a sohd phase pinned by the layer structure, is created. This low T phase 
has flnite helicity moduli and thus, zero resistance for any direction on the 
layer, while it has a nonzero vortex flow resistance guaranteed by zero helicity 
modulus in the perpendicular direction to the layers. Consequently, there is 
a transverse Meissner effect for a tilt perpendicular to the layers but not 
for any tilt parallel to the layers. With a slight change of the p-value, a 
structural transition possibly mediated by a unpinned solid occurs between 
different pinned solids and is reflected on the Tm{B) line in intermediated 
flelds as its oscillating 5-dependence [26]. This oscillating behavior has been 
suggested in YBCO data [55,56]. However, a description of p-dependences of 
such consecutive structure transitions is highly complicated and will not be 
given here. 

The isotropic form of the gradient terms in eq.(3.3) within the layers 
(in x-y plane) leads to a key insight on the physical picture in the liquid 
regime. It implies that, for high enough p values, the linear responses, such 
as the resistance, measured along the layers are independent of the relative 
direction between 5A (i.e., the current) and B. This is the essence of the 
so-called in-plane Lorentz force-free behavior observed [57,58] in tesla range 
of BSCCO where p > 1 is safely valid. Simultaneously, the in-plane isotropic 
form of gradients implies that, on cooling, the phase coherence lengths grow 
isotropically on the layers, while, as in the case B ± layers [34], the phase 
coherence perpendicular to the layers above is not sensitive to cooling 
and remains microscopic as a result of the (partial) LLL degeneracy. Since 
the phase correlation must be compatible with the positional correlation 
of vortices [25,26], the above-mentioned anisotropy appearing in the phase 
correlation must be also satisfied by the positional correlation. Hence, the 
observed in-plane Lorentz-force free behavior proves that, above T^, the po- 
sitional correlation first grows along the layers rather than across the layers, 
which is an opposite trend to an argument favoring a vortex smectic liq- 
uid [59]. Namely, the observation [57,58] is incompatible with assuming the 
intermediate phase [59]. 

It is not easy to describe in details the glass transition due to point 
disorder in this case B || layers. In ref.27, a high field behavior of glass 
transition line, corresponding to eq.(2.16) in B ± layers, was found and its 
interpolated behavior to lower fields was conjectured. 
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Applying a similar analysis to that leading to eq.(2.14) to the present 
case, the irreducible vertex /irr of xg — ~ lirr)"^, when A^cor ~ 0(1) 
^ 1, is found to have the form 

4r ^ ^exp(-2p), (3.5) 

where ji in the present case satisfies 

/x~exp(2/xo/4"'nr)) (3.6) 

which is identical with eq.(2.10), and the random flux A$ was again ne- 
glected. The e~^-dependence of eq.(3.5) is a reflection of /9a — 1 ~ e~^, 
i.e., of a weak spatial variation of and is also related to an exponentially 
small shear modulus of the pinning- free solid [60,26], while the exponential T- 
dependence in eq.(3.6) is a reflection of a 2D-like superconducting fluctuation 
weakened by a partial breaking, due to the layering, of the LLL degeneracy. 
Further, the p-insensitive melting line, eq.(3.4), is a consequence of eq.(3.6) 
which is also independent of p. 

Here we will assume the first order transition line Tm{B) to terminate 
at some p (> 1) and hence to have a upper critical point indicated as Pc in 
Fig.4. This is reasonable, because the above-mentioned exponentially small 
shear modulus in p > 1 likely results in a stronger random-pinning effect 
with increasing p. Then, even the first order transition in p < p^ may not be 
accompanied by an ordinary superconducting ordering signaled by /i — > 0, 
and eq.(3.6) becomes valid even below T^. Consequently, using the above 
expressions, one obtains a Josephson-vortex glass (JG) transition line 

^^^^^^ " MWT) V-^fd?-T^V-' A jj ' ('-"^ 

which has a similar T-dependence to eq.(2.26). Note that the prefactor de- 
creases with increasing 3D fiuctuation strength Of {(X VT) and that, due to 
the assumption Acor ~ 1, the A ^ limit cannot be taken in eq.(3.7). 
In the sketched phase diagram Fig.4, cq.(3.7) is useful just in p > pd. 
In pc < p < Pel, an exponential decay of A^cor resulting from a similar p- 
dependence of shear modulus will not be negligible with increasing p. Then, 
the prefactor of eq.(3.7) effectively diminishes in this field range, and the 
resulting i?jG(T), as described in Fig.4, becomes more flat. This phase di- 
agram Fig.4 will be compared with existing data in HTSC An a.c. suscep- 
tibility measurement for examining the onset of lock-in phenomena (i.e., of 
the transverse Meissner effect) in BSCCO has been performed [28] , and the 
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onset line of lock-in behavior found there has shown a relation similar to 
eq.(3.7) appearing in p > p^i in Fig.4. A detailed resistivity measurement 
has been performed in an optimally-doped BSCCO above 60 (K) [29] and 
has shown, at lower temperatures, a continuous vanishing of resistance along 
a flat curve consistent with Bjq{T) of Fig.4 but, near T^o, a discontinuous 
resistivity vanishing at Bjg{T) just on or below the disorder- free melting 
line Tjn{B). Further, recent resistivity data in 60K YBCO have suggested 
the presence of a remarkable slush regime and of a B\cp — 7 (T) above which 
Tm is roughly independent of p [56]. The occurrence of a Bi^p in p ~ 1 of the 
case B || layers is not surprising if the system is moderately dirty so that the 
VG transition line in p < 1 lies above the Tm(-B)-line. In fact, the "vertical" 
Tm{B) in p > 1 suggests that a situation with _Bicp and with the first order 
transition at higher p occurs more easily than in the case B ± layers. The 
continuous resistivity vanishing in 90K YBCO reported previously [61] might 
be a phenomenon below a -Bicp- 

In the case with disorder, the resistances for a current in all directions 
vanish simultaneously at Tjq{B). Further, the transverse Meissner effect 
signalled by a critical divergence of tilt modulus for a tilt across the layers 
is expected to occur for any p- value, i.e., even if the vortex lattice in the 
pinning — free case is a unpinned solid. A detailed analysis leading to these 
conclusions on response properties is seen in ref.27. 

4. SUMMARY AND DISCUSSION 

We have briefly explained the existing theory of phase diagram and phys- 
ical properties deep in the vortex liquid regime based on the LLL approxima- 
tion of the GL model. This theory can describe the vanishing behaviors of 
resistance in various situations in a way consistent with the phase diagram, 
while it is unclear to what extent the glass phases and possible transitions 
between them are described starting from the original GL model. 

In §2, a comparison of the theory with experimental data was done for 
YBCO. In BSCCO with stronger fluctuation and much larger anisotropy, 
it is known through the angular dependence [62] that the internal gauge 
fluctuation (leading to a magnetic screening changing features of interaction 
between the vortices) is no longer negligible in the low fields where the first 
order transition is realized, and hence, the present theory will not be directly 
applicable. 

Nevertheless, we note that, due to the strong dependence on the anisotropy 
of bI^^ in eq.(2.27), the absence of a lower critical point in BSCCO is not 
surprising. The position of the upper critical point and the presence of slush 
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regime in BSCCO are not clear yet to us. We simply guess here that, like 
in Fig. 4, the VG transition curve and its high field limit (corresponding to 
eq.(2.16)) in BSCCO should lie extremely below an extrapolated curve of 
Tjn{B) to higher fields. 

Recently, the presence of pinning-induced growth of \(Jxy\ near a glass 
transition was shown theoretically [63] and experimentally [64,65] and can- 
not be explained correctly based on the mean field vortex dynamics neglecting 
thermal fiuctuation [66]. This should be the case, because the present exten- 
sion of the fiuctuation theory [6] has explained the details of phase transition 
lines at which the linear dissipation vanishes. The glass fiuctuation contri- 
bution aG,xy to the Hall conductivity is expressed, as well as a^xx, by Fig.2, 
and it is expected that the magnitude of a-G,xy grows on approaching a glass 
transition with its sign, relative to that of (Xp^xy, dependent on the dimen- 
sionality of a dominant pinning disorder [63]. The phenomena in the case 
with only line disorder arc the best understood ones, and we find for this 
case that crG,xy ■ (XF,xy < [63,65] and that the ratio (Jcxy/ o'cxx (i-e., the Hall 
angle near the transition) does not vanish but seems to approach a constant 
at the transition [65,66]. Its details will be reported elsewhere. 
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c 



Figure 1: Diagram representing lu-r- The straight hne denotes the LLL prop- 
agator, the dashed hne is the "impurity" hne carrying A, and the semicircles 
imply vertex correction. 




(a) (b) 



Figure 2: Two typical Feynman diagrams contributing to aG,xx, in which 
the chain line denotes the next lowest Landau mode of ijj, and the hatched 
rectangle implies the correlation function (1.1). 




Figure 3: Schematic phase diagrams conjectured for B ± layers. The 
specific case (b) follows from the generic one (a) as a consequence of strong 
fluctuation. The sohd curve denotes the second order glass transition, while 
the chain curves include both the thermal first order line and the BrG melting 
line. 
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Figure 4: Schematic phase diagram for the case B || layers. Tm{B) and 
Tc{B) denote, respectively, the melting curve in the pinning-free case and 
the Hc2{T) crossover hue. The Bjq{T) curve and the first order transition 
line arc expressed by, respectively, the solid and chain curves. Any reflection 
of possible structural transitions below Tm{B) is not described here. 
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